New Mexico Tech (February 1, 2008) 



Noncommutative Einstein Equations 



Guglielmo Fucci and Ivan G. Avramidi 

New Mexico Institute of Mining and Technology 
Socorro, NM 87801, USA 
E-mail: iavramid@nmt.edu, gfucci@nmt.edu 



We study a noncommutative deformation of general relativity where 
the gravitational field is described by a matrix-valued symmetric two- 
tensor field. The equations of motion are derived in the framework of this 
new theory by varying a diffeomorphisms and gauge invariant action con- 
structed by using a matrix-valued scalar curvature. Interestingly the gen- 
uine noncommutative part of the dynamical equations is described only in 
terms of a particular tensor density that vanishes identically in the commu- 
tative limit. A noncommutative generalization of the energy-momentum 
tensor for the matter field is studied as well. 



1 Introduction 



The basic notions of General Relativity (GR) are closely related with the geomet- 
ric interpretation of second order elliptic partial differential operators describing 
the propagation of massless scalar fields without internal structure LI J- Let us con- 
sider a manifold M without boundary. Let / e C°°(M) be a smooth function on 
M and let L : C°°(M) — > C°°(M) be a second order partial differential operator of 
the form 

L = -a tlv (x)d fl d v + b u (x)d fl + c(x) , (1.1) 

with smooth real coefficients d", W and c. In particular, is a real smooth 
symmetric non-degenerate matrix at each point of the manifold M, that is, for any 
point x in M 

ct\x) = a v \x) , detcf v (x) ± . 

Since the matrix is symmetric, it has real eigenvalues and since the determinant 
is different from zero through the whole manifold M, it has constant signature. 
If the number of negative eigenvalues is 1, the operator L is called hyperbolic. 
Hyperbolic operators defined on M describe the propagation of massless scalar 
fields, in particular photons, over the manifold M. 

It can be easily shown that for the operator L to be covariant the matrix cf v 
should transform like a tensor of type (2, 0) under diffeomorphisms of the man- 
ifold M. This enables us to identify the matrix with a Riemannian metric 
guv over y[ Once a metric g^ v is defined one can construct all the geometrical 
quantities needed to develop the theory of General Relativity. In this theory the 
gravitational field is described by a symmetric non-degenerate tensor of type (0, 2) 
and its dynamics is described by Einstein's equations. 

A natural generalization of the GR theory leads us to consider the gravitational 
field described by a (NxN) matrix of metrics (g MV ) A B instead of a single metric (for 
details see [EDS)- Such a matrix-valued metric appear naturally if one considers 
the propagation of a field with an internal structure (a multiplet of fields rather 
than a single field). In this case the coefficients of the operator L, that is, the 
functions a MV (x), b M (x) and c{x), become matrix-valued. This follows from the 
idea that at short distances, or at high energies, the role played by the photon can 
be played by a multiplet of gauge fields Q]|. 

It is important to stress here that this approach for deforming general relativity 
is different from the ones proposed in the framework of noncommutative geom- 
etry (51|6l|7l|9l, where the coordinates do not commute and the standard product 
between functions is replaced by the Moyal product (S). 
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The purpose of this paper is to derive the equations of motion for the field a^ v 
that generalizes the role played by g^ v in the general theory of relativity. Since this 
model is a noncommutative extension of Einstein's general relativity we will call 
the corresponding equations of motions noncommutative Einstein's equations. 

The outline of the paper is as follows. First we develop the necessary formal- 
ism that will be used in the derivation of the dynamical equations of the theory. 
More precisely we will describe a generalization of various important geometri- 
cal object to endomorphism- valued quantities. Then we construct an action that 
is invariant under diffeomorphisms of the manifold M and gauge transformations 
and contains the matrix-valued scalar curvature. By varying, then, the action with 
respect to the new dynamical field a MV , we obtain the equations of motion. 

2 Matrix Geometry 

We will describe, in this section, the formalism needed in order to write the action 
for matrix general relativity following 01 |2j and derive the equations of motion 
of the theory. This formalism is a generalization of differential geometry, which 
is the natural language of general relativity, and it will be called Matrix Geome- 
try. The formalism that we are going to describe is related to the algebra-valued 
formulation of Mann and Wald flU. 

As already stated before, the main idea is to describe the gravitational field 
as a matrix-valued symmetric two-tensor field. Let M be an n-dimensional Rie- 
mannian manifold without boundary and let g MV be a metric tensor defined on the 
tangent bundle TM. It is well known that General Relativity is nothing but the 
dynamical theory of the metric 2-tensor field which is, basically, an isomorphism 
between tangent and cotangent bundles. In our model the metric 2-tensor field 
is replaced by a endomorphism-valued 2-tensor field a MV which represents an iso- 
morphism of more general bundles over the manifold M. The main idea here, sim- 
ilar to general relativity, is to develop a dynamical theory of this endomorphism- 
valued 2-tensor field cf v . This generality brings a much richer structure and con- 
tent to the model. 

Let V be an Af-dimensional Hermitian vector bundle over M, let ST = TM <g> V 
be the bundle constructed by taking the tensor product of the tangent bundle to 
the manifold M with the vector bundle V, and let Zf* = T*M (g> V, where T*M is 
the cotangent bundle to M. We consider, then, the bundle Iso(^, &*\ elements 
of which, are isomorphisms between 2T and In general the isomorphisms 
B : — > S7* can be identified with the sections of the bundle T*M <8> T*M <g> 
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End(V) and the isomorphisms A : SF* — > S? with the sections of the bundle 
TM® TM ®End(V). 

Let a be a symmetric self-adjoint element of TM <g> TM <g> End(V), that is, 

a"" = a^, {a 11 ")* = cf v . (2.2) 

This element is an isomorphism between the bundles 2?* and & if the following 
equation 

o">„ = iff , (2.3) 
with arbitrary if/ e ^7*, has a unique solution 

cf> v = . (2.4) 

This last requirement can be cast in another form; the element a is an isomorphism 
if the following equation 

cfbyp = b, v a vp = % ■ I , (2.5) 

has a unique solution b e T*M® T*M® End(V). 

There are some properties of the matrix b MV that need attention. The first prop- 
erty is the following: the matrix b^ v satisfies the following equation 

b% = b VM , (2.6) 

but it is not necessarily a self-adjoint matrix symmetric in its tensor indices. More- 
over, one can use a MV and b MV to lower and raise indices, although particular care 
is required in this operations because, in general, a MV and b MV do not commute and 
b M y is not symmetric in its tensorial indices d). 

We need, now, to develop some kind of invariant calculus. For this purpose, 
we introduce matrix-valued Christoffel symbols srf^ap that transform, under dif- 
feomorphisms of M, as connection coefficients, namely 



, d*" dx y dx 5 ^ , x dx' 11 d 2 x v _ 

St * a'B'(x ) = -st \ S (X) + • I , (2.7) 

PK ' dx v dx' a dx'P y V dx v dx' a dx'P 



where I is the identity matrix. Once we have a connection on M we can define a 
way to differentiate tensors. 

Let Tp(M) be the tensor bundle of type (p, q) on the manifold M. We define 
the following linear map 

£> : T« ® V — » Tf { ® V (2.8) 
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by 



p i 

J - J aH>v l ...v q u a l f J y 1 ,..v q t / i AaVvi,..v q nu ( r Vl ,., Vi _ l to M „.v q • > 

7=1 i=l 

This map is a well defined operator between T% <g> V and 7^ +1 <g> V. It is important 
to stress, at this point, that the linear map defined in (12.91 ) is not a covariant dif- 
ferentiation. In fact, it is easy to show that, because of the non-commutativity, the 
Leibnitz rule does not hold. However, in the commutative limit, the operator (12.81) 
reduces to an ordinary covariant derivative on the manifold M. 

In complete analogy with the ordinary Riemannian geometry, we shall define 
the matrix curvature and the matrix torsion tensors. Let (p e T*M® V, we compute 

D^Dyfa - DyDpfa = -R\ MV (f>A + T\ v T) A <p a , (2.10) 

where 

7?V = - d v $/ A aM + £f A ^ av - =g/V<^% > ( 2 - 1 1) 

and 

rV = ^V-^V (2.12) 
Once the matrix curvature (Riemann) tensor is defined we can construct the matrix 
Ricci tensor, namely 

K v = n a flav . (2.13) 

In order to write the action for matrix gravity, we need to introduce the matrix 
scalar curvature 'R. Since the metric (f v and the Ricci tensor are matrices, they 
do not commute in general and the definition of the scalar curvature, obtained by 
contracting the metric tensor with the Ricci tensor from the left, would be different 
if the contraction would be performed with the metric tensor on the right. In order 
to avoid this choice, we use a symmetrized definition of the matrix-valued scalar 
curvature as follows 

n= l -(a^ v n llv + n^ v ) . (2.14) 

We need, now, to relate the connection coefficients s^ p iiy to the metric ten- 
sor a /Jv . We impose a compatibility condition similar to the one in Riemannian 
geometry as 

+ ^V*^ + = • ( 2 - 15 ) 

The above equation has the following solution in a closed form 0]| 

= ^b A(T {a ay d y a p<T - a^dyd™ - crd y a ap + S ap(T + S p<Ta + S crpa )b m , (2.16) 
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where S is an arbitrary matrix-valued tensor that satisfies the symmetry property 



S a„v = _ r vp _ (2 _ 17) 

The matrix-valued tensor S is related, in the general case, with the torsion T in 
(I2.12I ). Moreover, in the commutative limit, the tensor S reduces exactly to the 
torsion. It is important to notice that in matrix geometry the connection (12.161) is 
not symmetric in the two lower indices even if the tensor S vanishes. In the rest 
of the paper we will assume, without loss of generality, that the tensor S vanishes. 

In order to write an action for the model under consideration we need a gener- 
alization of the concept of measure. As a guiding principle, any generalization of 
the measure /u has to lead, in the commutative limit, to the ordinary Riemannian 
measure detg^ v |. Moreover the measure /i is a density depending only on the 
metric a MV and not on its derivatives and transforming, under diffeomorphisms of 
M, as 

lu'(x') = J(x)/u(x) , (2.18) 

where 

7 W = det(^). (2.19) 

A definition of the measure /i which is a straightforward generalization of the 
Riemannian measure, is the following [OQ 

fi = ^Tr v p, (2.20) 

where p is a matrix- valued scalar density, which can be defined, for example, as 
follows 



r % 



exp(-a^ v ) . (2.21) 

Then p only depends on the metric a and transforms in the correct way under 
diffeomorphisms of M. 



3 Variation of the Action 

We construct the action functional for the field following [1J. This functional 
has to be invariant under both diffeomorphisms of M and gauge transformations. 
The infinitesimal form of these transformations is 

V v = [w,aH, (3.22) 
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and 

5^ v = ^cT , (3.23) 
where a> is an element of the algebra of the gauge group 

8JJ = Co , 

and £ is the generator of the infinitesimal coordinate transformation 

<s f jc" = ^ = -e . 

We can construct an action functional for field a MV that satisfies the properties 
described above, by using the matrix-valued scalar curvature, defined in (12.141) . 
and the measure (12.201 ). A good candidate for the action is the following 

S M GR(a) = -^— f dx2-Tr v (pK). (3.24) 
!6nG J M N 

It is worth noticing that because of the cyclic property of the trace, the relative po- 
sition of p and the scalar curvature is irrelevant, moreover it is easily shown that 
the action functional (13.241) is invariant under the diffeomorphisms (13.231) and un- 
der the gauge transformations (|3.22l) . Since the action has the invariant properties 
discussed above, the currents associated with the symmetries satisfy the identities 

M^h^J^ - <3 - 25) 

and 

6S 



*•! °- (3 ' 26) 
where the (13.251) is the current generated by the invariance with respect to dif- 
feomorphisms and (13.261) is the current generated by the internal (gauge) symme- 
try. The above identities represent an endomorphism- valued generalization of the 
Noether identities, which are related, in the usual theory, to the contracted Bianchi 
identities. 

Now that we have the action functional we can derive the equations of motion 
for the field a^, which is the main goal and the main result of the present paper. 
These equations will be matrix-valued and they will constitute a generalization 
of the ordinary Einstein's equations that we will call noncommutative Einstein 
equations. In order to find the dynamics of the model we vary the action (13.241) 
with respect to the field a MV considered as independent variable, namely 

a ^ _> a ny + SaT . 
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By doing so we obtain, for the variation of the action, the following 

1 f 1 

6S = S(d" + 6a MV ) -S(a MV ) = dx-Tr v (3uvSa MV ) , 

\6nG Jm N 



(3.27) 



where Q MV is some matrix valued symmetric tensor density. Then, of course, the 
desired equations of motion are 



Guv = . 



(3.28) 



It is important to notice that the matrix-valued tensor density (|3.28l) has to coincide 
with the Einstein tensor in the commutative limit, more precisely we need that, in 
the commutative limit, the following relation holds 



1 



l 



N 



(3.29) 



Our main task, then, is to find the explicit form of the equations of motion 
that result from the variation of the action (13.241) . In all the calculations that will 
follow the order of the terms is important, unless explicitly stated, due to the 
matrix nature of them. 

First of all, we rewrite the action in a more explicit form which is more suitable 
for the subsequent variation, namely 



MGR. 1 



(a) 



i r i 

— dx— Tri 
!6nG J M N 



(3.30) 



By varying the terms in (13.301) with respect to the independent field a MV , and by 
using the cyclic property of the trace we get 



1 C 1 

dS MGR (a) = —— dx—Tr v 
\6nG J M N 



(3.31) 



where the curly brackets { , } denote anti-commutation, namely {A, B} = AB + BA. 

From the expressions (12.1 II ) and (12.131 ), we can evaluate the variation of the 
matrix- valued Ricci tensor, more precisely we have 



8<R 



From now on, for simplicity of notation, we set 

b^ v = {p,cr\. 



(3.32) 



(3.33) 
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By substituting (13.321) in (13.311) , and by using the cyclic property of the trace we 
obtain 



SS mgr(«) 



i r l 

— — dx— Tr, 
nG J M N 



\6nG 
1 

2 



Sp<R + ^W MV ,p}Sa» v + i^V? (1 ( ( U/" /n .) + 



1 



By integrating by parts and by collecting similar terms we get 
£Smgr(«) 



(3.34) 



f dx^Tx v \5 P n+U% v ,p}8a^- l -(B^\ a -B^^ A aA + 
IbnG j M N L 2 2\ 



We can rewrite the last expression in a more compact form, namely 

1 . .. ... 1 . _ 1 



.(3.35) 



55 M gr(«) = TT~^ \ dx Tr Tr ' 
16^-G J M /V 



(3.36) 

where the matrix- valued tensor densities C MV a and LP have the explicit expression 



C»\ = {a^,p, a - p^\ A \ - p[aT, £t\A ~ [p, ^ pa W + 



- {p, W v Aa , a^] } - a^W v M ,p\ + 2{p, a^W\ a 



and 



D 1 = {a">, v - stf p vp p) - [p, ^ pv ]a pv - [p, 



(3.37) 



(3.38) 



It is worth noticing that in the commutative limit, or, in other words, when all the 
matrices commute, the tensor densities O lv a and D M are identically zero, and the 
variation of the action 55mgr simply reduces to the standard result of the general 
theory of relativity. 

We can write, now, the variation of the connection coefficients. By using the 
expression (12.161 ), and by noticing that 



Sb^ = -b lip (6(f a ')b (rv , 



'MP 
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we obtain the following 



£/ a Av 6a v % u + -b A(T 6a ay (d y a p<T )b pu + 



- -btoStfridytf^bpp - -bxadcHidycDb^ + -b A(T a ay {d y 5a p<T )b pp + 



1 

T 
1 

2 l 



-b^idySanbp, - -b A<T a^(d y 6a pa )b p , 



(3.39) 



Once we have the explicit expression for the variation of the connection coeffi- 
cients, we can evaluate the last two terms that appear in the variation of the action 
(13.361) . We start with the first of the two 



- u- 

L Jm 

- kS' 

z Jm 



rfjcTr v (C"V^°™) 



dx Tr 



^ a ts v C pv a b MA Sa^ + b 0v C pv a ^\ A 6a A/3 + 



- -{d y a p<T )b pv C pv a b p(T 5a ay + -{d y a a<T )b pv C pv .b^ + 



1 

2 [ 



+ -(dyCf ^bpvC^ab^da^ - -b pv C> lv a b li(r a (rY {d y 8a f " T ) + 
+ ^b pv C p \b p(T a py (d y 6a a<T ) + l -b pv C p \b p<T a^{d y 5a pa ) 



(3.40) 



where in this last expression we used the cyclic property of the trace. 

We introduce the following definition, which will be useful in order to simplify 
the notation, 

i> P = bfrfTJow . (3.41) 
By using the above definition, the expression in (13.401) can be rewritten as follows 

- \ [ dxTr v (C pv a 6^ a MV ) = 



dx Tr , 



M 



^% y c? p F paA 8a A l 3 + F pap a py ^ a yA 6a^ + 



-(d y a (Tp )F pa(T 6a ay + ^(d y a a<T )F pa(T 6a py + ^(d y a'"')F riil ,6a" y + 
1 1 1 



1 

2 [ 



-F patT a ay (d y 6a p(T ) + — F pa(r cf r (d y 6a a " r ) + -F pa<T cT(d y 8cr) 



, (3.42) 
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where the first two terms in the last expression has been derived by using the 
relation 

^V^aV 5 ** = ^V^'M^V 5 *** = sd a Py a yp F paX Sa^ . (3.43) 

By integrating by parts and by relabelling dummy indices we find the final expres- 
sion for (13.421) . namely 



z Jm 



dx Tr viC^adsrf 



V) "of' 
L Jm 



dx Tr i 



a riyP F 4- F nPV ( 
PyU r p a y t r p a pli s 

6a A " + 



yA 



+ 



+ -id. 



F p A<r F Acrp F * <TpA 



Ay 



y a p<T . (3.44) 



For the last term in the variation of the action (|3.36l) . we use similar arguments 
which lead us to the expression (13.441) . In this case we introduce the following 
definition: 

Gpap = bfkiybv . (3.45) 
By using the definition above and the cyclic property of the trace we obtain 

\ f dxTr v (irS^% a ) = 
z Jm 

- IL 

L Jm 



dx Tr i 



/a P A8a^ + 



-(d y a p(T )G pa(T 6a or + -(dX^G^Sa^ + -(d^^G^Sa^ + 



- -G pa(T a ay (.d y 6a pcr ) + -G pacT a py (d y 6a acr ) + -G pa<T cf y (d y Sa pa ) 
By integrating by parts and relabelling dummy indices we get 



. (3.46) 



z Jm 



dx Tr yip's*?* m ) 



z Jm 



dx Tr i 



a py G ya x + Gp ay a 



py rxf a 



pA 



+ 



-jidf}a P<T ){G p Aa- - Gxcrp - Ga-pxj 



- - d 



GpXcr — Gxa-p — Ga-pX \d 



6a A/3 + 
Sa p(T . 



(3.47) 
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It is worth noticing that in the above expressions, (13.441) and (13.471 ), the ten- 
sor densities F and G always appear in the same combination. This observation 
justifies the following definitions 



X, 



pAa 



F pAcr F Xcrp F (T p A 



(3.48) 



and 



IpAcr 



GpAcr — Gxcrn ~ G, 



Aap 



T crpA 



(3.49) 



By using the two definitions above we can rewrite the arguments of the traces in 
(13.441 ) and in (13.471) respectively as 



r^-^ 1 a^S^ pv 



^ a p y a w F paA + F pap a py ^% x - -{dptf^X^ + 



+ -d 7 (X^) 



da*? , 



(3.50) 



and 



ir 

"2 



'V" 



cFGyrt + G^^ a pA - -{d^ p<T )Y pA(T + 



1 



da"' 3 



(3.51) 



By combining the results (13.501) and (13.511) we obtain the expression for the last 
two terms in the variation of the action, namely 



- ^C"WV + \^^ a m = \{^ a p P a py {F yaA - G yaA ) 



+ 



HFpay - GpryWsTpy + ^aniYpAcr ~ X pA(r ) + 



(Y Ap p - X Ap fj)a 



■py 



6a^ 



(3.52) 



4 Noncommutative Einstein Equations 



With the expression (|3.52l) for the last two terms in (|3.36l) . the variation of the 
action has the form (13.271 ) which is suitable for the derivation of the dynamical 
equations of the model. Before writing the complete dynamical equations, we 
will simplify further the expression (13.521) . 
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The definition (13.481) gives a linear relation between the matrix- valued tensor 
density X and a particular combination of matrix- valued tensor density F, a similar 
linear relation between Y and G is given in (13.491) . By using simple tensor algebra, 
it can be easily shown that those relations can be inverted, namely we can write 

FpAcr - -^(XAcrp + X apA ) , (4.53) 

and 

G pA(T = -jiYzrp - Y<t P a) ■ (4.54) 

By substituting the equations (14.531) and (14.541) in the expression (13.521) we obtain 
the following 

-ic"V*V + \D p 6^% a = ^% p a py [(Y aAy -X aAy ) + (Y Aya -X Aya )] + 
+ i(Y ayA — X ayA ) + (Y y p a — X y p a )]a py £/ a pA + 
+ (d p an(Y pA a - X P Mr) - d y [(Y Apf3 - X^a^Sa* 3 . (4.55) 

We can see, in the last formula, that the tensor densities X and Y enter always 
in the same combination. It is useful, therefore, to define the following tensor 
density 

H pvp — Y pV p Xfiy P . (4.56) 
With this last definition we can rewrite (14.551 ) as 



icv^V + \w5^ a pa = i 



'% p a^H aAy + ^ a ^R Xya + 



8a Ap . (4.57) 



+H ayP a py ^ a pA + H yfia cP zf* pA + (d fi a pcr )H pAxr - d y (H Apl3 a py ) 

By using the compatibility condition (12.151) , we can write that 

d p a p(T = -£/ p y/3 a y<T - ^yp-aT , (4.58) 

moreover we obtain that 

- d y (H Ap/3 a py ) = -{d y H ApP )a py + H^^^cH + K^sT^ . (4.59) 
Since H pvp is a tensor density, we can write 

DyH Ap p = d y H Ap p - s$ a ' x r H ap p- - srf a ' py H Aa p - £/ a /3 y H Apa - .s^ a ya FL A pp . (4.60) 



12 



By using the results obtained in (14.581) , (14.591) and (14.601 ) we can express (14.571 ) as 
follows 



-(D y H App .)a p ? - W a Ap ,H ay p]a py - H ay ^ a Ap ,a py ] - [^ a py ,H Aa ^ + 



-W a ya ,H ApP \a p T - ^% y [H Apa ,a p ^6a^ .(4.61) 
At this point we introduce the operator P defined as 

PyH Ap p = DH Ap p + [£/ a Ay , H ap p\ + \sf* py, H Aa p\ + [£/ a p y ,H Apa ] + [^/ a ya ,H Apl s\ . 

(4.62) 

By using the last definition in (14.611) one obtains 

- \o\5^% v + i/XWV = -^ a WA a^H aAy + 2H aAp a^^\ A] + 
-{P y H Ap p)a^ + W% y , H Apa ]a^ - H ayP W a Ap , a??] - ^% y [H Apa , a"] W" . 



(4.63) 

We finally have all the ingredients that we need in order to write the dynamical 
equations of the theory. Now we only have to find an expression for the variation 
dp. The definition of p is given in (12.211) . and its variation can be straightforwardly 
evaluated as follows 

6p = - f % f ds e^-Wfir^-MO , (4.64) 

J 7T2 Jo 

R" 

where 

A(0 = a^ v . (4.65) 

Once we have the expression (14.641) for the variation, we can use the cyclic prop- 
erty of the trace to write that 



Tr v (6p R) = Tr , 



- f d 4 f 

J 7T2 Jo 



' ds e- sA ®<Re-V- sW %& 



6a^ v . (4.66) 



By combining (14.661) . (14.631) and (13.361) we obtain the noncommutative Ein- 
stein equations in absence of matter, namely 

= , (4.67) 
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where 

Q, v = \{p, Kp,} + T pv + \^\ p] a p m avy + l -H aAv a p ^ a W] - l -(P y H ppv )a p ^ + 
+ l -W a vy , H^cf* - l -H ayv Y^% P , a py ] - ^\[H upa , a p ?] , (4.68) 
is the noncommutative Einstein tensor, is defined by 

T, v = -[ d 4C ds e-* A ®<Re-^%{ v , (4.69) 

J 7H Jo 

R" 

and the tensor density H has the explicit form 
H aAy = b^8\D»-C»\)bw-b x MyD»-C»\)b m -b yv (8\D»-C»\)b lxA . (4.70) 

These equations are the main result of the present paper. One can show that the 
first two terms in the equations (14.681 ) represent a straightforward generalization of 
Einstein's equation to endomorphism- valued objects and the rest of the terms can 
be considered as a genuine noncommutative part which is not present in Einstein's 
equation. It is interesting to note that the pure noncommutative part is completely 
described by the tensor density H Mvp defined in (14.70I ). 

Moreover the equation (14.671 ) satisfies the requirement (13.29K which, in words, 
expresses the necessity that our model reduces, in the commutative limit, to the 
standard theory of general relativity. In fact, the trace of the pure noncommutative 
terms vanishes, because of the presence of the commutators, and the first two 
terms just give 

iTr y |i{p,^ v } + rJ= Vstav-^v*l • (4.71) 

For an arbitrary matrix algebra the equation (14.671) becomes more complicated 
than the ordinary Einstein's equation due to presence of the new tensor density 
H^p. We mention, now, a particular case in which (|4.67l) simplifies. The formal- 
ism used so far deals with geometric quantities which are endomorphism-valued, 
namely they take values in End(V). By choosing a basis in the vector space V we 
can represent End(V) by means of matrices. Let us suppose that the algebra under 
consideration is Abelian, in this case all the elements commute with each other 
and the tensor density H^ vp vanishes identically and the equation (14.671) becomes 

Kv ~ 2 = • (4-72) 
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Therefore, in case of a commutative matrix algebra, the equation of motion of our 
model have the same form as Einstein's equation, with the only difference that 
(14.721) is matrix-valued. 



5 The Action for the Matter Field 

In order to have a complete theory for the gravitational field we need to describe 
the dynamics of the matter field in the framework of matrix general relativity. The 
main idea is to extend the general results of classical field theory. We will consider, 
in the following, the dynamics of a multiplet of free scalar fields propagating on 
a manifold M. We can construct an invariant action by using the matrix valued 
metric and the measure p. A typical action is 

SmatterO, <f) = j f dx {- (dtf, [p, a^}d v (p) - (if, [p, Q)(f)} , (5.73) 
^ Jm 

where ( , ) denotes the fiber inner product on the vector bundle V, and Q is a 
constant mass matrix determining the masses of the scalar fields. The equations 
of motion of the scalar fields are then obviously 

[-d lt {p,ar}dv + {p,Q)]<P = 0. (5.74) 
The complete action of the gravity and matter is described then by 

S (a, <f) = S M Gr(<3) + Smatter(a, <p) . (5.75) 

By varying the above action with respect to one obtains the noncommutative 
Einstein equation in presence of matter 

g, v = SnGNT^ , (5.76) 

where T MV is the matrix energy-momentum tensor defined by 

letter 

By using the explicit lagrangian (15.731) for the matter field, we obtain the expres- 
sion for the energy-momentum tensor 

T pv = - [{p, d M ip (8) dM + AV + N MV ] + (ju^v), (5.78) 
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where the explicit form of M^ v and N^ v is obtained by using the variation of the 
scalar density p in (14.641 ), namely 

AV = " f^r f ds e- sA ^a a ^d a <p®dw}e-V- s W%Z v , (5.79) 

J 7T2 Jo 

R" 

and 

= " f ^ f ^ ^{fi, <P 9 <p}e- {1 -^%{ v . (5.80) 

J 7T2 Jo 

R" 

It is worth remarking, here, that the above formula (15.781) for the energy- 
momentum tensor 7^ v reduces, in the commutative limit, to the standard result, 

e.g. m. 

6 Conclusions 

The main idea of this new model is to describe the gravitational field by a multi- 
plet of gauge fields with some internal structure. For this purpose the metric field 
g flv , which describes gravity in general relativity, is replaced by a matrix-valued 
2-tensor field a^ v . This allows the model to have a much richer content in de- 
scribing gravitational phenomena. A more general geometric picture is developed 
by allowing the metric to be matrix-valued. Most of the geometric quantities, 
used in describing gravity, can be generalized to be endomorphism- valued. In this 
framework it is possible to introduce an action for the gravitational field which 
is diffeomorphisms and gauge invariant, that leads, after performing the variation 
with respect to a^, to the modified (noncommutative) Einstein equation. It is in- 
teresting that the noncommutative part of the modified equations only depends on 
a specific tensor density H^p and on a linear combination of its commutators. 

With an explicit expression for the noncommutative Einstein equation, it will 
be possible to study some particular simple solutions of (|4.67l) : for example a 
static and spherically symmetric solution, which would describe the gravitational 
field outside a spherically symmetric massive body, or a spherically symmetric ho- 
mogeneous solution, that would describe noncommutative cosmological models. 
These simple examples are very interesting and we plan to study them systemati- 
cally elsewhere. 

We would like to make a final remark. In our model all the geometric quan- 
tities that we need to develop the formalism are endomorphism-valued. Once a 
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basis for the vector bundle V has been fixed, we can represent elements of End( V) 
by matrices. Of course the description of physical phenomena has to be inde- 
pendent from the particular realization of the representation. This is, ultimately, 
related to the gauge invariance of the theory. We believe that by an opportune 
choice of gauge, namely an opportune representation of End(V) by matrices, the 
dynamical equation (14.671 ) could be simplified further. The search for such partic- 
ular gauge, if it exists, requires further studies in matrix differential geometry and 
matrix general relativity. 
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